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We present calculations ol the charge density profile, subband occupancy and ellipsometry spectra 
of the electron gas at the LaA10;i/SrTi03 interface. The calculations employ self-consistent Hartree 
and random phase approximations, a tight binding parametrization of the band structure and a 
model for the optical phonon of SrTiOs. The dependence of the spatial structure and occupancy of 
subbands on the magnitude of the polarization charge at the interface and the dielectric function is 
determined. The interface-confined subbands may be labelled by the symmetry (xy, xz, or yz) of 
the Ti d-orbitals from which they mainly derive. The xy-derived band nearest the interface contains 
the major proportion of the electronic charge, but a large number of more distant, slightly occupied 
xy-derived bands are also found. Depending on the magnitude of the polarization charge, zero, 
one, or two xz/yz derived subbands are found. When present, these xz/yz bands give the dominant 
contribution to the long-distance tail of the interface charge. The response to applied ac electric 
fields polarized parallel and perpendicular to the interface is calculated and the results are presented 
in terms of ellipsometry angles. Two features are found: a dip in the spectrum near the LO feature 
of the STO phonon and a peak at the higher energy. We show that the form and magnitude of the 
dip is related to the Drude response of carriers moving in the plane of the interface while the peak 
arises from the plasmon excitation of the xz and yz electrons. The relation of the features of the 
subband occupancies and the in-plane conductivities is given. 



PACS numbers: 71.27.+a, 73.40.-c, 73.50.Mx, 78.30.-j 

I. INTRODUCTION 

Transition metal oxide materials have been investi- 
gated extensively because of their remarkable electronic 
properties, including high Curie temperature magnetism, 
metal insulator transitions, spin and orbital order and 
superconductivity^]. In bulk materials these phases can 
be controlled by doping, pressure and chemical compo- 
sition but the degree to which material properties can 
be varied is limited. Progress over the last decade in 
pulsed laser deposition and molecular bean epitaxy has 
provided new ways to control the properties of novel ma- 
terials, by enabling the formation of atomic precision su- 
perlattices and interfaces unobtainable with bulk crys- 
tal growth methods 3 -"- Artificially created oxide sys- 
tems have been reported to exhibit a range of strik- 
ing behavior inc ludin g metal-insulator transitions^^, 
superconductivity^ 3 -^, anc j magnetisrrP^NASI The very 
short length scales and low volume carrier densities 
in these systems mean that characterization-of crystal 
structure, of defects, of electronic structure and of car- 
rier density- is a key issue. 

Perhaps the most widely studied class of artifi- 
cally created strongly correlated systems is the "ox- 
ide interface" J 3 -^ where in appropriate circumstances the 
junction between two insulating materials can support a 
two dimensional electron g as wit h interesting p ropert ies 
including superconductivity^^] anc j magnetis m 1 15 * 17 * ^. 

The most widely studied system is the interface between 
LaA10 3 (LAO) and SrTi0 3 (STO^P but variants in 
which the LAO is replaced by LaTiOsEM] GdTiOaEH 
or vacuurrpS r the STO by KTaO^ have also been 
studied. A wide variety of experiments and theoretical 
calculations suggest that the electron gas is confined to 



within ~ 10 nm or less of the interfac e 1 14 * 25 * ^, resides in 
the STO in subbands derived from the Ti £ 2g symmetry 
d-orbitals^ 3 25 26 29 ^ and is induced by the polar catas- 
trophe scenaricP^SlEfl] However, the details of the value 
and spatial extent of the charge density as well as the 
subband structure and occupancy are not clearly under- 
stood. 

In the LAO/STO system the carrier density inferred 
from Hall effect and other transport measurements is typ- 
ically one order of magnitude less than the number theo- 
retically expected from polar catastrophe^, although high 
energy photoemisson measurements yield numbers more 
nearly consistent with theory^. Other contributions to 
the charge density besid es the polar catastrophe mecha- 
nism have been propose d 33 * 34 ! and a key role for defects 
has been documented^. 

Inferences from the mag netic field dependence of the 
superconducting propertie d 7 * 14 * indicate that the spa- 
tial extent is ~ 10 — 20nm whereas band theory 
calculations^^ place most of the electrons substantially 
closer to the interface. Finally, the t2 g orbitals into 
which the electrons are believed to go have a substan- 
tial anisotropy in their dispersion, with one of the three 
bands dispersing much less strongly in the direction per- 
pendicular to the interface than the others. The relative 
occupancy of the different orbitals therefore plays a cru- 
cial role in the spatial structure of the electron gas. 

In this paper, we employ a minimal tight binding 
model to describe the subband structure and occupancy 
of electrons at a polar interface. Our model includes 
charge self-consistency at the Hartree (self-consistent 
screening) level and the spatial and frequency structure 
of the dielectric constant of the host material as well as 
the orbital degeneracy of the ti g levels that host the rel- 
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evant electrons. The subband structure and occupancies 
are computed as a function of the magnitude of the polar- 
ization charge at the interface. We find that most of the 
added charge goes into bands derived from d xy orbitals 
which lie very close to the interface, but that an impor- 
tant contribution to the density of charges far from the 
interface comes from electrons in the lowest-lying d xz / yz - 
derived bands and that the form of the long distance 
charge distribution reflects the spatial structure of the 
lowest subband wave function. The dielectric tensor is 
then calculated within the random phase approximation 
(RPA) and is used in a solution of the relevant Maxwell 
equations to obtain expressions for the response to in- 
cident radiation polarized parallel and perpendicular to 
the interface. The results are expressed in terms of re- 
flectivity spectra and ellipsometry angles. The relation 
of the ellipsometry angles to the charge density, subband 
occupancies, and in-plane conductivity is determined. 

The rest of this paper is organized as follows. Section[IT| 
introduces the tight binding/Hartree model and presents 
the theoretical formalism needed for the calculation of 
the dielectric tensor and reflectivity. Section |TTlj presents 
the subband structure, charge density, and orbital oc- 
cupancy following from the solution of our equations in 
different physically relevant cases. Section [IV| is devoted 
to the analysis of the ellipsometry angles. Section [V] is a 
summary and conclusion. 



II. FORMALISM 

A. Tight-binding model and self-consistent Hartree 
approximation 

In this subsection, we introduce a tight-binding model 
describing the electronic structure of electrons at an ox- 
ide interface. While our specific motivation is the elec- 
tronic properties of the LAO/STO interface, the basic 
physics we address is that of electrons in ^-derived 
bands of a transition metal oxide material confined by a 
Coulomb potential arising from a charge sheet. We there- 
fore believe our results apply also to the STO/vacuurrPSl, 
GTO/STO-^and (with appropriate modifications) to the 
LAO /KTaCP^l systems. Extension to systems with e g de- 
rived bands would be straightforward but is not consid- 
ered here. 

Bulk LAO and STO are insulators with band gaps of 
about 5.6 eV and 3.2 eV, respectivel y ^ 36 1 37 1 . They crys- 
tallize in forms which are small distortions of the stan- 
dard ABO3 perovskite structure and we use the simple 
perovskite notation with axes (001), (010), (100) aligned 
along the B—O bonds to indicate directions. We consider 
the most common experimental situation, in which an 
(001) interface of a very large STO crystal is capped with 
some number M of layers of LaA103 in the same crystal- 
lographic orientation. When viewed along the (001) di- 
rection, LAO is polar and if the interface is such that the 
Ti02 plane of the STO adjoins the LaO plane of the LAO 



then basic electrostatic consider at ions^IEI] suggest that 
in the idealized situation of perfectly aligned defect-free 
materials, the electric fields caused by the p olar cha racter 
of the LAO induce an interface electron ga d 3 * 21 * 35 4 

The areal charge density of the induced electron gas 
vanishes if the number of capping layers M is less than a 
critical value of order 2 — 4 and for M greater than this 
critical value the areal charge density increases, eventu- 
ally asymptotic to th e value 0.5 electrons per in-plane 
unit cell. TheorjE^H and experiment!^ indicate that the 
critical M at which charge first appears, the thickness 
beyond which the induced charge becomes close to the 
theoretical asymptotic value and indeed the asymptotic 
value itself depend on the dielectric properties of host 
(STO) and overlayer (LAO) material, on the chemistry 
of the free surface of the LAO (which controls the en- 
ergetics of the compensating charges) and on density of 
charged defects such as O vacancies or interstitials. For 
our purposes it is enough to model the polar interface 
as a semi- infinite block of STO, terminated in the (001) 
direction by a sheet of positive areal charge density. The 
magnitude en /d 2 of the areal charge density (d 2 is the 
area of a unit cell.) is a parameter of the model, and we 
consider a range of values between uq — and no = 0.5. 

In the polar catastrophe scenario, for thick enough po- 
lar over layers, a charge of no is transferred from the STO 
region to the outer boundary of the polar overlayer. The 
compensating electronic charge density — n is trapped 
near the interface by a potential which has a contribution 
V — Aire 2 riQz/ed 2 from the charge at the outer bound- 
ary of the polar over layer. This potential is reduced 
by the scale dependent dielectric parameter of STO and 
by screening from the interface electrons. In this situa- 
tion the in-plane crystal momentum is a good quantum 
number while an orbital-dependent subband structure 
describes the motion perpendicular to the plane. We 
treat this situation using a self-consistent Hartree ap- 
proximation. The band alignments of LAO and STO 
are such that it is energetically prefe rable to a ccommo- 
date the trapped electrons in the STd 3 | 26 | 27 J3 ^, in bands 
derived from the t2 g {d xy ^ yz%xz ) orbitals of the Ti. We 
model these bands in terms of a tight binding model 
with three orbitals per unit cell and only nearest neighbor 
hopping. The nearest neighbor-only hopping assumption 
implies that the orbitals are not mixed on the single- 
particle level, and moreover implies that electrons in the 
xy-derived bands move only in the xy plane etc. This 
drastically simplifies the ensuing theoretical treatment. 
In the actual band structure, small farther neighbor hop- 
ping terms exist which mix the orbitals. However, the 
effects of these terms are small and do not affect our 
main conclusions, which have to do with charge distri- 
butions and average optical responses, in any important 
way. 

With these considerations, we write the effective 
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Hamiltonian of the system as 

ff(k) = E ^(k) c&ck 

k 

?= {xy,yz,xz} 

+ E (^W + ^jc^ct , (i) 

k i 

where k is the in-plane momentum, i and j are plane in- 
dices, I denotes orbital indices, *L(k) represents the hop- 
ping between Ti atoms (orbitally diagonal, within our as- 
sumptions) , V/ 1 [n] is the self-consistent Hartree potential 
arising from the electron density n, and represents 
the external potential from the polarization charge. 

Because we are interested in low electron densities per 
orbital we approximate the hopping matrices t\j (k) as 

i*?(k) = td 2 |k|% 

#/(k) = td^Sij - *(<%_i - 2S iS + 5 ij+1 ) , 
tf?(k) = td 2 k 2 J l3 - t(5ij-i - 2% + , (2) 

Here we have made the small k approximation to the 
dispersion and observe that corrections to the orbital- 
diagonal nearest-neighbor-only hopping approximation 
appear only to higher order in k. d is the distance be- 
tween two Ti atoms (any difference between the in-plane 
and transverse Ti-Ti distance is irrelevant for our pur- 
poses and is neglected). We set t = 0.34 eV. 

The spatial dependence of the potential arising from 
the polar charge and the self-consistent field of the elec- 
trons is complicated by the large and strongly scale- 
dependent behavior of the dielectric function of ST(j 39 * 40 4 
Within our tight binding approach we associate electrons 
with Ti sites and the dielectric function becomes a link 
variable, depending on the distance of a given Ti — Ti 
bond from the interface. The potential arising from the 
polarization charge then may be written as 

K^^nog— , (3) 

with tjj+i the dielectric constant appropriate to the re- 
gion between the Ti-site in plane j and the Ti site in 
plane j + 1. The values of the dielectric function 

are considered as parameters that arc fitted by compar- 
ing the resulting density with existing DFT results^. 
Similarly, the Hartree potential Vf 1 ' is written as 

where rij represents the charge per unit cell in the j th 
Ti0 2 plane, given as rij = (£, < =cy ) = (£ qi C'^C 1 ^}. 
Charge neutrality implies no = n i- We note that the 
Hartree potential has no dependence on in-plane momen- 
tum in this approximation. . 



The charge densities are obtained from the solution 
of the self-consistent Hartree equations; this is greatly 
simplified by the orbital-diagonal and nearest-neighbor- 
only nature of the Hamiltonian. The basic equations are 

E [*-,( k ) + <%(Vf [n] + V^)} ^-(k) = £ L^(k) , 

3 

(5) 

for V'ij(k) and e l kn , the wavefunction of the nth eigen- 
state of orbital I and the corresponding eigenvalue, and 

nt= E M4JIV>L(q)l 2 (6) 

nq 

I — {xy,yz,xz} 

with the Fermi-Dirac distribution function, iV(£q„) and 
chemical potential chosen so that the total electronic 
charge per in-plane unit cell equals the polar charge no- 
Three features of the equations warrant comment. 
First, the spectrum involves both bound (to the inter- 
face) and unbound eigenfunctions. Second, because in 
our formalism the xy orbital disperses only in the plane 
parallel to the interface the n th xy band involves carriers 
in plane n (measured from the interface) and has min- 
imum energy s^L 0n = V^ at + V^. At long distances 
from the interface the binding energy tends exponen- 
tially to zero with Thomas-Fermi decay length set by 
the long-distance behavior of the dielectric constant and 
by the hopping t. The very large value of the STO dielec- 
tric constant means that this length is in practice very 
long. Third, the xz/yz eigenfunctions are delocalized 
over many planes, and as we shall see in reasonable cases 
only one bound-state subband is occupied. 

B. Current correlation function and dielectric 
tensor 

In this subsection, we calculate the current correla- 
tion function using self-consistent linear response (RPA) 
theory applied to the Hamiltonian of the previous sec- 
tion. We relate the obtained response functions to the 
transverse dielectric tensor and use them to calculate the 
reflectivity of the electron gas at the interface. In de- 
riving the current correlation function, we assume that 
the impurity scattering may be treated within the relax- 
ation time approximation, i.e. may be modeled by a fre- 
quency and momentum-independent scattering time 
The resulting formula gives valid results both in diffusive 
(cj <C 1/t) and collision-less (w 3> 1/r) regimes. 

The response functions depend on the polarization of 
the incident light, in other words on the direction of the 
applied electric field with respect to the plane of the in- 
terface. For in-plane electric fields, screening effects are 
negligible as are interband transition matrix elements 
(the small further-neighbor hopping terms neglected in 
our theory would give weak and unimportant interband 
terms). The response is therefore a standard Drude re- 
sponse in which each occupied subband contributes in 
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parallel. Because the yz (xz) band does not disperse in 
the x (y) direction, the x-direction (y- direction) current 
has contributions only from the xy and xz (xy and yz) or- 
bitals. Therefore for the in-plane response we can expect 
that we can obtain information concerning the orbitally 
resolved electron density per in-plane unit cell, n 



xy 



or n 



xy 



"yz, 



as well as the in-plane conductivities. 



The electronic response from out-of-plane electric fields 
arises only from the xz/yz bands because in our approx- 
imation the xy orbitals have only in-plane dispersion. 
Unlike the case of in-plane electric field, interband tran- 
sitions (principally to unbound states) are important. In 
the physically relevant situation of an effectively infinitely 
thick STO substrate, the unbound states form a contin- 
uum. In our actual calculations we use a finite system 
with a large number N (typically 40-80) of STO layers 
but the modest impurity-induced broadening that we in- 
clude means that our calculated results are a continuum 
whose form has negligible dependence on N and broad- 
ening. 

Computations of the optical response require the gauge 
invariant current operator, which is obtained by inserting 
the vector potential A into the Hamiltonian in the usual 
minimal coupling way and expanding to zeroth and first 
order in A leading to "paramagnetic" and "diamagnetic" 
terms which are expressed in second quantized notation 
as 

2" 



Jp(q) = 2etd J2 2(p-f)J^ r 



, C p-qm /C P m ' 



pmm' 

J^(q) = 2etd J2 2<&U 



p-qm"-"Pm ! 



pmm 



r ||,dlo 



q) 



E eL,C^J& m A!(q, W ) 

kmm' 



e 2 t 



+^m>ii>mi+l)Al , i(q) 



(7) 



where the factor 2 in each equation is from the spin de- 
generacy, c qm is an annihilation operator for m th sub- 
band eigenstate with momentum q satisfying C qi = 



the A is the vector potential, c is the 



speed of light, the || and _L denote the in- and out of- 
plane directions and the superscripts p and dia represent 
the paramagnetic and diamagnetic contributions to the 
current, respectively. Because the hopping and hence the 
current is orbital-diagonal we have not explicitly denoted 
the orbital index here; this information is included in the 
eigenvalue index m runs from 1 to 3iV with N the num- 
ber of Ti02 layers included in the calculation including 
the orbital index. The current in z direction is properly 
defined as a link variable; we use i+\to denote the link 
between the i th and (i + l) st plane. The current depends 
upon the wave functions; this information is encoded in 
the functions 



(8) 



and 



i+ i 1 

Qm'm = TTli^Pm'i+l^mi — ^mH^Pmi+l) 



2d 



ltd 



/ j m r 



(9) 



To obtain the expression for Q we have used Eq. ([5| . An 
important simplification following from our assumption 

of nearest-neighbor hopping is that both F % m , rn and Q m i 2 m 
are independent of in-plane momentum q since eigenvec- 
tors tpmi for eigenvalue e qm are the same for all value of 
q and the difference between two eigenvalues with the 
same q is independent of q. 

We adopt the gauge in which the scalar potential van- 
ishes and the vector potential is related to the electric 
field as A = -^-E and initially consider the clean limit 
in which the impurity scattering time r — ¥ oo so that 
the problem is translation-invariant in the plane of the 
interface. Linear response theory then yields the para- 
magnetic and diamagnetic contributions to the gauge- 
invariant conductivity tensor as 



1 3 



Aie 2 t 2 d 



2D , , W - 
kmm' 

(k+^)F j ,+ 

\ 2 mm 



' £m'k ^mk-fq ^~ IT) 

- o , 



* mm 



±.dia 



m£xy 



2ie 2 t 
cod 

^ t n mF mm XX 
m£xz 

1 \- 



E n ™ F m m yy 



m£xz 
or yz 



,(io) 



where the factor 2 in a - ?- is due to the spin degeneracy, 
N2D is number of unit cells in the plane, rj is an in- 
finitesimally small positive number, and n m is defined as 
( 2 S q c qm c qm) including spin degeneracy. 

We note that in the limit of q — > 0, the off-diagonal 
tensor components such as xy vanish since <7?-(q, uS) — 
— ofj{— q, Ui). Moreover for xy bands, the clean limit 
assumption which will be relaxed later means that the 
conductivity comes only from diamagnetic term since 



F': 



In this limit the conductivity tensor is 
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reduced to 

, Aie 2 t 2 d tt-^ 

m<Exz k 



2ZY mrno f p (£fcmo ) 
~2 



A 2 

mmo 



2^ ^mm a ^m m + + 4Q m m Qm rnZZ 

Aie 2 t 2 dn, 



IT) 



(A:™ ) 2 

o F mmQ Fl nam {xx + yy) + 2Q m m Q mo mZZ 



2ie 2 i 
I6ie 2 t 2 dn 



^ ] ^rnm -^m Q m{^ : ^ : VV) 



m^Lxz 



A, 



-ZZ , 



(11) 



where we use the symmetry between xz and yz orbital, 
take the index m to run over the eigenstates of xz orbital 
and assume that only the lowest yz and xz band, denoted 
as mo, is occupied (this will be discussed in more detail 
below). Therefore n mo becomes the number of electrons 
in yz (or xz) ban ds per in-plane unit cell (n mo = n yz = 
n xz ). In Eq. <[TTJ) we define nf = Emexy n ™ F L m whicn 
is the number of electrons per in-plane unit cell in 
band located at i th layer and 



2A r . 



j2 - A mm + "7 



(12) 



and have used the identities 

S-F l F i F j 

^ij memo / j mmo mom ? 
m 

M^moi+lVW + ^m iVW+l) 

= 8W 2 



^cmmo ^mom 



(13) 



exploiting = ip m i. Combining paramagnetic and dia- 
magnetic terms, we obtain the conductivity in the limit 
q — » as 



2ie 2 i 



2+ r 



<%nf (:H + yy) 



+n mo £ ( ^f°^T„ M + 1 ) F: nmg Fl om (xx + yy) 

m£xz 

+n mQ 



A mmo (w 2 - A 2 imo + it?) ' 



(14) 



We note that the resulting conductivity can be separated 
into "local" (i.e. layer-diagonal) (superscript L) and non- 
local (layer changing) (superscript NL) contributions as 



a-ij = 






2ieH 


[(n 


ujd 






±JfL 
ij 


2ie 2 t 


"2 


iod 


3 ( 



which are given as 



+ n* z )xx + (nf + nf )yy)} 



c u 1 



'•m a P mo ( w )^rr 



j F m m (,XX - 



yy) 



+ 



n 2 ± 



8n mo tui 



l d 2 Q 



,(^ 2 



A 2 

mmo 



+ if})' 



Z£ F p 



■into 3 



"'■0 



4?rf 

where we define nf 



O 2 i i 

"-L, ,pm 7 l +2 7J+2 55 
. . U/J mo rnra ^ ra ^ ra 



(15) 



F l 

mo^ fflflmo 



plasma frequency Si'' 2 = 87re 2 £(nf 
cal conductivity, and introduce f2j_ = 



and in-plane 

+ n,y z )/d for lo- 
8ire 2 tn mo /d and 



Zmmo — - Ej=i -^mmo for non-local conductivity with 
the use of Eq. ^ . For electric field independent of layer 
index (as is the case for radiation polarized parallel to the 
interface), the first term of is zero (^2- F^ mo = 0) 
and thus for optical transitions, we can ignore the first 
term. Then the structure of the conductivity tensor be- 
comes much simpler: the in-plane response is given by 
the sum of the Drude conductivities of each layer and 
the out-of-plane response is from inter-subband transi- 
tions. The oscillator strengths in the conductivities are 
characterized by the plasma frequencies ilj' and Q±. 

We now include impurity scattering within the re- 
laxation time approximation. 41 From the relation = 



e°lS ir 



stant without interface electrons, we obtain the trans- 
verse dielectric tensor with relaxation time 1/r as 



eij(u) = e (ui)Sijl 



) 2 r 

— — ——dijixx + yy) 
ui{ui + i/t) 



If 



i+i 



m£xz 



where 



L m {u) = 



2tA, 



A 2 

mmo 



- iuj/r 



(16) 



(17) 



In the regime where u) S> 1/r, the dielectric response in 
the direction perpendicular to the plane is characterized 
by a series of inter-band peaks broadened by Suj ~ l/2r. 
On the other hand, for uj <C 1/t, we have diffusion poles 
where the imaginary part of dielectric function goes to 
zero in the limit of to — > 0. In Eq. (16 1 we find that 



dielectric tensor in the direction perpendicular to the in- 
terface is finite in the limit of to — > 0, while for in-plane 
direction it is proportional to \/to. This comes from the 
our assumption that electrons at the interface are con- 
fined in the direction perpendicular to the interface but 
may propagate freely in the plane of the interface. 
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The calculations of the reflectivity spectra in the next 
section require the inverse of ef?. From the expression 



in the Eq. ( 16 ), we can write the inverse of the dielectric 
tensor as 



O 2 \ T 7 1+ 2 7 
|_ / ^m^mmo^'. 



'2 r/J ^ 2 

mom 



_l_r)4 \ 1 t 7*+ 2 7*1 + 2 r 7*1 + 2 7J+2 



^"^-L ^ ^ $ mm 1 -^mi $m\7n' -^m' 



O 2 <? , f , 



m^m' 



mi^xz 



Sij + ni 



mmo 



L- 1 -n 2 ± s 



z 



mom' ' 

(18) 



where we define 



5,, 



LI 



e^ 2 



We introduce the Coulomb matrix 

J rnrn 



limit of g — » is related with S, 

,. 27re 2 
hm 



' as: 



(19) 

(g), which in the 



lim ?; - - lim \ p -9«*l»-il f-J 



47re^ 



5,, 



With the definition e°[Ae 2 



(20) 



kj 



(5y we have 



where 



n, 



= 2n r . 



mm'£xz 



tL^ 1 - 2n r . 



(21) 



(22) 



and Eq is a bare charging energy defined as Aire 2 /d. In 
the following section we show that (e°(cj)) 2 [Ae Z2: (a;)]~- 1 
functions as an effective dielectric constant in the expres- 
sions for the reflectivity. We also note that the poles of 
density correlation function H mm ' in Eq. (21) represent 
the inter-band plasmon excitations. These are shifted 
from the inter-band transition energy due to the Coulomb 
interaction. 

If the perturbing electric field is very smooth compared 
with the system size, it is reasonable to approximate 
[Ae*»]r/ ~ EylA^Hly 1 = «Ae-»», implying 

(e°) 2 «Ae- 1 )) =E° C £ (m \z\m)U mm> (m'\z\mo)(23) 

mm' ^xz 



usin S E^-nmo = Hi iF mm = { m o\ z \ m ) which is the 
dipole matrix element between the ground and excited 
states. Thus in the long-wavelength limit, we can ex- 
pect that the structure of (e°(w)) 2 ((AeJ 2 , 1 (a;))) depends 
on plasmon poles and the dipole matrix element between 
occupied and excited states. 



n = 





1 --- 


no = 
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5 — 




FIG. 1: (color online) Imaginary part of the (e ) 2 ((AeJ z 1 (w)}} 
for polar charge density of 0.1, 0.3, and 0.5 electrons per in- 
plane unit cell. For simplicity we set as a constant e° = 
lim^^oo £g(w) = 5.1 where £%(uj) is the dielectric constant of 
bulk STO. The scattering rate is set to 1/r = 40 meV for all 
three cases. 



In order to investigate the characteristics of plasmon 
excitation, we calculate (e°) 2 ((AeJ 2 1 )) from the solution 
of self-consistent Hartree equation (Eq. (|5|) by assuming 
typical values for the static dielectric constant {ejj+i} 
which will be discussed in Section |III| For simplicity we 
set e°(ui) as a constant value (e° = 5.1) that is the high 
frequency limit of the dielectric constant of STO. Fig. [I] 
shows the imaginary part of the (e°) 2 ((AeJ 2 1 )) for polar 
charge density no of 0.1, 0.3, and 0.5 electron per unit 
cell. We can see a well-defined peak for all three cases. 
The peak energy becomes higher and the peak width be- 
come broader as polar charge density increases. In fol- 
lowing sections, we will show that this plasmon excitation 
produces high energy features in ellipsometry angles. 



C. Reflectance and ellipsometry spectra 

In this subsection, we obtain formulae for reflectiv- 
ity by solving Maxw ell's equation using electromagnetic 
Green function d 42 * 43 ! and the dielectric function obtained 
in the previous subsection. We consider two types of 
incident radiation: s-polarized, where the incoming and 
outgoing electric fields are both parallel to the interface 
plane and the orthogonal p-polarization in which we have 
both in-plane and out-of-plane components of electric 
fields. 

From the Maxwell equation, we can write the differ- 
ential equation for the electric field and current density 
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V x (V x E) = 



18 

~cdi 



1 <9E 

c ~dt 



(24) 



By introducing the dielectric tensor e{uj,r, r') = 
e (u), z)S(z — z') + Ae(r, r') and transforming to frequency 
space, we can rewrite the equation as 



-V x V x +^-e (cj,z) 



E(w,r) 



/" dVAe(w,r,r')E(w,r') 



(25) 



In Eq. (25 1, the contribution to the dielectric tensor 
of the interface electrons is denoted as Ae(w,r,r / ). It 
is in general non-local in space but vanishes in regions 
where the electron density vanishes. We assume that 
the other contributions, such as phonon modes in STO, 
are local in space and are included in e° and define the 
bare dielectric constant of LAO/STO heterostructures as 
e°(oj,z) = e° L (oj)l for z < (i.e. in the LAO region) and 
e°(u>,z) = e° s (cj)i for z > (i.e. in the STO region). 
In addition, we can write the electric field components 
as E(w,r) = E(cj, z)e zC *' p ~' lut with in-plane wave vector 
Q and position p where |Q| = -^/e^wsin^/c with inci- 
dent angle 9. In the rest of this section we will omit the 
frequency dependence in the expressions for the electric 
fields and dielectric constants for simplicity. Solving the 
Eq. ( 25 1 is similar to solving a scattering problem where 



we have non-local dielectric function in the place of the 
scattering potential and we can write the general solution 
asP 



E(z) = E°(z) 



d Zl dz 2 G(z,Zi)Ae(zi,z 2 )E(z 2 ) , (26) 



where 



-V x V x + — e°(z) 



and E°(z) satisfies 



-V x V x 



-e\z) 



E°(z) =0 



(28) 



The explicit expression of G and E° for s- and p-polarized 
lights are presented in appendix A. 

For s-polarized light, the incoming and outgoing elec- 
tric fields are both perpendicular to the plane of incidence 
and we take the direction of the electric field to be the y 
direction. Eq. (26 1 then becomes 



E V (Z) 



dzidz 2 G yy (z, zi)Ae y y(zi,z 2 )Ey(z 2 ) 



(29) 



Given the system size L ~ 10nm, the wave number times 
system size kL is estimated as \kL\ ~ |47r\/e°ai?§^T§| <C 
1, where otp is the fine structure constant. This means 
that neither the Green function nor the electric field 
changes significantly in the region where intcrfacial elec- 
trons reside. Thus, we can approximate the Eq. (29) as 



E y {z) = E° y (z) - ^G yy (z,0)((Ae y y))E y (0) , (30) 
where 

«Ae w » = Jdz ldz2 Ae yy ^z 2 ) 



uj(uj + i/r) 



(31) 



and flu is defined as J2i ^f 2 or 8ire 2 t(n xy + n yz )/d. We 
note that in the clean limit only diamagnetic terms con- 
tributes to the approximated dielectric tensor since the 
non-local terms give a contribution of the order of (kL) 2 . 
Therefore, we expect a Drude response of the interface 
electrons for s-polarized light. By substituting z — on 
both sides of Eq. (31 ), we can obtain the E y (0) and write 



the electric field as 



E y (z) 



e lqz + r°e- tqz 



cuMl + rj) ((Ae w ))(l + r°) 



iqz 



e lqz + r®e- tqz | 1 



2iqd 

Zo J) 



((Ae w )) 

(l+r°)qd 



(32) 



where r° is the reflectivity without interfacial electrons 



G(z, z') = lS(z - z') , (27) which is 



^-tt , and o and k are wave vectors in z-direction 
for z < and z > defined as ojy/e^ cos9/c and 
u) yCg — \ \ sin 0/c, respectively, for an incident angle 9. 
In Eq. (32), the relation (1 + r° s ) 2 = 4g 2 r° is 
used. From the definition of the reflectivity E y (z) — 

iqz 



iqz we have 



2iqd 



((Ae OT » 



; s 1 



(l+rg)qd 
2ie°, cos 2 9 



(33) 



We note that the reflectivity becomes —1 in the low fre- 
quency limit and r° in the high frequency limit as ex- 
pected. 

For p-polarized incident light, we have both in-plane 
and out-of-plane components of electric fields. Assuming 
that the electric field is in the plane of incidence (xz 



plane), we can write^ 
E x (z) = E° x (z) 



dzi dz 2 [G xx ( z , zi ) Ae xx ( z\ , z 2 ) E x ( z 2 ) 



+ G xz (z,z 1 )Ae zz (z 1 ,z 2 )E z (z 2 )] 
E z (z) = E° z (z) 

, ,2 



dzidz 2 [G zx {z, z 1 )Ae xx (z 1 , z 2 )E x (z 2 ) 



+G' zz (z, zi)Ae zz (zi, z 2 )E z {z 2 )] 
-jj|-y J dz x Ae zz {z, zi)E z (zi) . 



(34) 



Using that the displacement field, defined as D z (z) = 
e°(z)E z (z), is continuous across the interface, we can 
simplify Eq. |29|) in terms of the displacement field as 



E x {z) 



E° x (z) 

, ,2 



dz\dz 2 



G xx (z,Zi)Ae xx (zi,z 2 )E x (z 2 ) 



-G xz (z, zi)e°(z 1 )Ae zz {zi 1 z 2 )D z (z 2 ) 
D z (z) = D° z (z) - 6 (*l U ( d Zl dz 2 



x G zx (z, z\ )Ae xx (z 1 ,z 2 )E x (z 2 ) 

-G' zz {z, z 1 )e°(z 1 )Ae-^(z 1 ,z 2 )D z (z 2 ) 
where we have used 

e°{z)Ae zz (z,z 1 ) = - / dz 2 Ae zz (z, z 2 )[e~ 1 ] zz (z 2 , z x ) . 



, (35) 



(36) 

By applying the same approximation as in the s- 
polarization case (Eqs. (30 1 and (31 )), we obtain the cou- 
pled differential equations 

E x (z) = E x '(z)-~G xx (z,0)((Ae xx ))E x (0) 

+ ^G xz (z,0)e s ((Ae z - 1 )) J D z (0) 



D z (z) = D° z (z) 



e° L G zx (z,0)((Ae xx ))E x (0) 



+ U ^e° L G zz (z,0)e° s ((^: z 1 ))D z (0) , (37) 
in which the dielectric properties of the materials are en- 



coded in the two functions ((Ae xx )) (defined in Eq. (31 )) 



which captures the in-pla ne P rude conductivity and 
((Ae zz )) (defined in Eq. ( 23 ) ) which represents the 
plasma excitation. Substituting z = leads to coupled 
equations for E x (0) and D z (0) that can be easily solved. 
By substituting the solutions back into Eq. (37), we can 
get the reflectivity r p defined as 



E x {z) 



r v e 



-%qz 



(38) 



Then the reflectivity expressed in terms of averaged di- 
electric tensor ((Ae xx )} and ((Ae zz )) is 

n ^((Ae xx )) + ^^ {{ Ae^) 



q -%^((Ae xx )) 



^°(l+rg 
2iq 



1 + 



liqd e" L 
f u _ f u f u „2_ f Q1 



l + ^l((Ae xx )) 



fc 2 ((A e ^)) + Q 2 ( e °) 2 ((A e - 1 ))) 

Q*de° L (l+r°) 



2 i<i 



'((Ae^ 1 )) 



(39) 



where r p = J q+e o k and we have used the identity (1 — 

r° P ) 2 = {(1 + r°)(eW° s 9)} 2 = r ^UA^ From 
r s and r p , we obtain the ellipsometry angles and A 
defined as 



and 



* = tan 1 (\r p /r s \) 
A = arg(r p /r s ) . 



(40) 



(41) 



III. RESULTS: ELECTRONIC STRUCTURE 

In this section, we present the electronic structure ob- 
tained from the solution of the self-consistent Hartree 
equations, with a focus on the dependence of the sub- 
band structure and the charge distribution on the value 
of the polar charge no- 

A solution to the self consistent equations requires val- 
ues for the static dielectric function between the dif- 
ferent planes. We estimate these by fitting our calcu- 
lated electron density to that obtained from DFT cal- 
culations including full lattice relaxation^. We obtain 
{e ii+1 } = {75, 350, 1000, 2000, 3000, 5000, 5000, ...} where 
we assign 5000 for i greater than 6 because our results 
turn out to be insensitive to the precise values of the di- 
electric constant, especially far from the interface, as long 
as these are large. The increase of the dielectric constant 
with increasing distance from the interface may be under- 
stood as a nonlinear effect arising from the electric field 
dependence of <PP: near the interface the internal electric 
fields are large, thus reducing e whereas at longer dis- 
tances the fields rapidly decrease, allowing e to approach 
the very large value appropriate to bulk STO. 

The calculations of Ref. [29] were performed for a situa- 
tion corresponding in our notation to polar charge equal 
to 0.5 per in-plane unit cell. For smaller polar charges one 
expects weaker electric fields and therefore larger e. In all 
of our calculations, however, we used the same effective 
dielectric constant for all polar charge densities since we 
find that the change in the dielectric constants fitted for 
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smaller polar charge density is not drastic and the char- 
acteristic features of electronic structure do not crucially 
depend on the detail of dielectric constants. Moreover, 
fixing e allows us to isolate the effect of changing polar 
charge on the electronic structure. 



(a) 



V/ 



n = 0.1 



0.02 
(b) 



0.01 



n = 0.1 




10 20 30 40 50 60 

z/d 

FIG. 2: (color online) Left panels: Electronic structure of 
interface electrons for the polar charge of (a) 0.1, (c) 0.3, 
and (e) 0.5 electrons for unit cell. Inset: magnification of 
lower lying yz and xz bands. Right panels: Charge density 
profiles corresponding to the three situations shown in the left 
panels. Inset: Fermi surfaces for the three cases. The amount 
of electrons occupied in xy and xz/yz bands are about 0.07 
and 0.03, respectively, for the polar charge of 0.1 electrons 
per unit cell, 0.17 and 0.13 for 0.3 polar charge, and 0.24 and 
0.26 for 0.5 polar change. 



Fig. [2] shows the band structure and charge density 
distribution of interface electrons for three different po- 
lar discontinuities. The self-consistent Hartree equation 
is solved for 80 Ti02 layers and we found that as long as 
the number of layers was taken to be greater than 40 the 
results are essentially independent of layer number. (As 
the number of layers goes to infinity negligibly small oc- 



cupancies of higher xz/yz subband may appear.) For the 
no = 0.5 case the results are consistent with published 
band structure calculations^!. 

The right hand panels of Fig. [2] show the total and or- 
bitally resolved charge densities. As expected, the charge 
is confined to the near interface region and the occupa- 
tions of the xy and xz/yz orbitals are inequivalent be- 
cause of the symmetry breaking arising from the inter- 
face. There are three spatial regimes of charge distribu- 
tion: a very large peak (arising physically from the xy 
orbitals) coming from the plane nearest the interface, an 
intermediate regime (z < 30<i) with a slow spatial decay 
and an appreciable contribution from the xz/yz orbitals 
(whose charge density in the very near interface regime 
is suppressed by wave function effects), and a far regime 
with a very small charge density, exponentially decaying 
with z and indistinguishable from in the Fig. [2] The 
bulk of the charge is accommodated within ~ 5 — 10nm 
(z ~ 10 — 25d) of the interface, in agreement with exper- 
iments. The spatial extent of the charge distribution de- 
pends on the magnitude of the polar charge, being largest 
for smallest polar charge and reflecting the strength of 
the confining potential. For all of the cases we have ex- 
amined, about 30% percent of the electrons ( 50% of the 
xy electrons) are located on the first layer. 

The left hand panels display the subband energies as 
functions of in-plane momentum in the x and y direc- 
tions. The d xz / yz bands are distinguished by their lack 
of dispersion along y and x respectively. For polar charge 
less than n ~ 0.04 (not shown) only xy bands are occu- 
pied. For larger no, only one xz and one yz band is occu- 
pied except for a narrow range of no very near 0.5 where 
there is a very small occupancy of the second xz and 
yz subband. Thus the z dependence of the xz/yz con- 
tribution to the charge density shown in the right hand 
panels arises almost entirely from the z dependence of the 
subband wave function. On the other hand a very large 
number of xy bands are very slightly occupied, and it is 
these bands which give rise to the very small Thomas- 
Fermi tail of charge density extending very far from the 
interface. 



IV. ELLIPSOMETRY ANGLES 

In this section, the optical response of the interface 
electrons is presented in terms of ellipsometry angles, 
which are a convenient experimental quantity because 
they permit a straightforward subtraction of the sub- 
strate and overlayer effects. In particular, even in the 
absence of an electron gas the optical properties have 
features associated with reflection from the LAO-STO in- 
terface. This structure is complicated by the very strong 
optical phonon of STO, which leads to strong frequency 
dependence of the STO dielectric constant, which is such 
that at one particular frequency the optical constants of 
LAO and STO nearly match and transmission through 
the interface becomes nearly perfect. Consideration of 
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changes in ellipsometry angles (defined more precisely 
below) allows these effects to be efficiently normalized 
out. The characteristic features of ellipsometry angles 
are identified and we will relate them with physical prop- 
erties of the system such as charge density, orbital occu- 
pancy, and scattering rates. 

In section II C, we derived expressions for the reflectiv- 
ity of s- and p- polarized incident light. For s-polarized 
light, the electric field is purely transverse and in-plane, 
so in the experimentally relevant long- wavelength limit, 
r s represents the in-plane Drude response of the inter- 
face electron gas. For p-polarized incident light with 
both in-plane and out-of-plane components of electric 
field; the breaking of translation invariance in the out- 
of-plane direction means that charge fluctuations are in- 
duced and the longitudinal response is relevant: the re- 
flectivity contains both Drude and inter-band plasmonic 
responses. Even in the absence of interface electrons r a 
and r p have structure associated with the STO optical 
phonon. In the absence of interface electrons, forming the 
ratio r p /r s removes this structure leaving a smooth be- 
havior. Thus, by comparing the ratio |r p /r s | with the in- 
terface electrons to the one without electrons, we can dif- 
ferentiate the responses of interface electrons from those 
from optical phonon. Ellipsometry experiments deter- 
mine ^> = tan _1 (|r p /r s |) and A = arg(r p /r s ) simultane- 
ously. Thus, calculating the difference A 'J between the 
ellipsometry angles with and without interface electrons 
enables us to investigate the properties of that electron 
gas at the interface. 

Fig. [3] presents our results as difference in ellipsometry 
angles A 'J with and without interface electrons for three 
different polar discontinuity of 0.1, 0.3, and 0.5 electrons 
per unit cell. There are two main features in the ellip- 
sometry spectra. One is a dip (highlighted in the left 
inset of Fig. [3] (a)) approximately at the plasma edge fre- 
quency ujp e where the real part of dielectric constant of 
STO becomes the same as the dielectric constant of LAO 
(assumed as the same as vacuum dielectric constant) and 
the other is a relatively broad peak at (highlighted in the 
right inset of Fig. [3] (a)) high energy related to the elec- 
tronic plasma excitation. The position of the high energy 
peak is seen to be roughly proportional to the number of 
electrons in the yz (xz) bands for the two larger charge 
densities. For the polar discontinuity of 0.1 the peak po- 
sition is affected by proximity to the STO LO phonon 
feature. A detailed analysis of the peak position will be 
given later in the section. 

In order to understand the relation between ellipsom- 
etry angles and the physical properties of the system, 
we consider that the electronic contribution to the re- 
flectivity is small compared with unity and thus the de- 
nominator of Eqs. (33) and (39) are approximated as 
one. The approximation is valid in broad frequency 
ranges satisfying both qd({Ae yz / xz }) / cos 2 8 <C 1 and 
qd((Ae~z)) / cos 2 <C 1 since the thickness of the con- 
ducting sheet is small relative to the wavelength of light. 
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FIG. 3: (color online) (a) Ellipsometry angles for the polar 
charge of 0.1, 0.3, and 0.5 electrons per unit cell with incident 
angle of 75°. Inset (left): magnification near the plasma edge 
where Re[e?,(w)] = Re[eg(w)]. Vertical line indicates bare 
interface plasma edge frequency to pe . Inset (right): magnifi- 
cation of the higher energy peaks, (b) Real and imaginary 
part of the dielectric tensor contributed by interface electrons 
for no = 0.3 



In this limit we have 



iqd(l +r°) 5 
2e° L cos 2 9 
iqd(l - r°) 2 



2e£ 



«Aey}) , 



«Aen» 



e°r sin 2 1 



A - e° sin 2 



; (e^) 2 ((A e - 1 )) 



(42) 



(43) 



where we define ((Ae\\)) = ((Ae xx )) = ((Ae yy }} using the 
cubic symmetry. 



and r p we 



With the approximated expressions for r 
consider two simplifying limits. First, we consider that 
the electronic contributions to the reflection relative to 
that from the mismatch of the two bulk dielectric con- 
stants is small and only take into account the effect of 
interface electrons up to the linear order. Then the ratio 
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r p /r s becomes 



P iqd 



(1 



„0\2 



r°(l 



el sin 2 i 



e° T sin 2 



(r°) 2 cos 2 6> y 
r(e^) 2 «Ae- 1 }) 



((Ae,,)) 



(44) 



the approximation used in Eqs. (42 1 and (43) is valid 
except very near grazing incidence (8 ~ 90°). More- 
over, we find that the in-plane dielectric tensor dominates 
especially near the STO plasma edge w. 



0.11 eV, 

which is also common for no = 0.1 and hq = 0.5 cases 
(not shown). Near uj pe , both the real and imaginary 



part of e% is small (Re[e|g] 



l,Im[e°] 



0.1) 



We expect that at high frequencies (w » 1/r) ((Aey)), 
essentially, Drude conductivity becomes small whereas 



and in this case the density correlation function II is 
roughly proportional to the inverse of the Coulomb ma- 
trix. Then the out-of-plane dielectric tensor can be writ- 



(e° s ) 2 ((Ae zz )) has a peak at an appropriate plasma fre- ten as -qd(e° s ) 2 ((Ae zz 1 )) ~ qL which is much smaller 
quency. In the low frequency limit (u> <C V T )) the con- 
tribution from ((Aey)) dominates due to the 1/uj factor 
in ((Aey)) whereas (e^) 2 ((Ae zz )) is finite as uj — > 0. 



In the opposite limit, which may occur at the 



than one. Thus, we expect that the dip at tu pe is due 
mostly to the in-plane conductivity. Since the real part 
of the bare dielectric constant of STO and LAO coincides 
at m ve and with its small imaginary part we can expand 



at which Re[e°] = Re[ejg], we have vanishing r° and Eqs. (42) and (43) up to the linear order of Im[e^] and 



if Im[ejg] is small. If the contribution from the in- electronic dielectric tensor 



terface electrons satisfies Im[eg] <C qd((Aeii)) / cos 2 9, 



qd(e%) 2 ((Ae- z 
written as 



z )) / cos 2 9, up to linear order, r p /r s can be 



JIm[e|]_ 
Ae° L cos 2 9 



M((Aeu)) 
2e° cos 2 9 



^( e °) 2 ((Ae- 1 ))/((Aey)) 



i(2cos 2 9- l)Im[e°] iqd((Aen)) 



ie° L cos 2 8 



2ei 



2ie° r cos 2 9 . n 
\r? 



r° cos 2 i 



gd«A C |,)) L p 
Vr s co S 2 9(e°) 2 ((Ae- z 1 ))/((Ae ll ))} 



(46) 



(47) 



(45) 



where we ignore the out-of-plane dielectric tensor and 
have used the bare reflectivity at co pe 



where the contribution from the out-of-plane excitation 
is expressed as a ratio of out of plane to in-plane dielectric 
tensor. 

In order to analyze the feature in the ellipsometry an- 
gles, we first examine the dielectric tensor ((Aey (w pe ))) 
and (e° s ) 2 ((Ae~ z (u pe ))) multiplied by qd/ cos9. Fig. 
(b) shows the real and imaginary part of ((Aey)) an 
(eg) 2 ((Ae zz )) for no = 0.3. First of all, we can see that 
overall magnitude of dielectric tensor is well below 1 so 



jlm[e|]_ 
4e^ cos 2 9 



i{2cos 2 9- l)Im[e£] 
Ae° L cos 2 9 



(48) 
(49) 



Since Im[e^] is in the same order of qd((Aey)), Eqs. (44) 



and (|45j) are not applicable. However, we can get some 
insight by evaluating the magnitude of r p /r s 



(2 cos 2 (9- l)Im[e°] 



2g pe c/0|| cos 



2q pe dQ^ cos 2 9/t 



1/2 



Im[e°] 



(50) 



where q pe is defined as the wavevector at uj p 



for small fl 2 and change its sign around 



For 9 < 45°, we can show that d\r p /r s \/ dVl 2 > 
regardless of the value of ft 2 and thus we expect that 
l r p/ r s| — at Lo pe is always positive and monotoni- 



cally increase as the density n xy 



increases. On the 



other hand, for 9 > 45°, d\r p /r s \/dfl 2 becomes negative 



(l-2cos 2 0)« + l/r 2 )Im[e^] 
2q pe dcos 2 9 



(51) 



assuming uj pe 3> 1/r. Thus, as we increase the polar dis- 
continuity, A 1 ^ will increase monotonically for 9 < 45° 
and for 9 > 45°, it will decrease until f2| ~ f2 2 c . We can 
also evaluate the effect of varying the scattering rate 1/r. 
In the regime 1/r <C uj pe , increasing the scattering rate is 
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FIG. 4: 

near oj v 



(color online) The structure of ellipsometry angles 
. (a) Magnitude of the dip at the plasma edge fre- 
quency as a function of n xy + n yz . The dots are from the 
simulated ellipsometry spe ctra and the solid line is obtained 
from the Eqs. (331 and (39 1 while ignoring the ((AeJ z }}. The 
dashed line is obtained by Eq. ( |50[ ). (b) Comparison of el- 
lipsometry angles with the different scattering rates for polar 
charge of 0.3 electrons per unit cell around the plasma edge. 



equivalent to changing fi 2 — ► fijj(l — l/(w pe r) 2 ). There- 
fore, the magnitude of A\I> will decrease as the scattering 
rate increases as long as 1/r < ui pe . Moreover, Eq. (50) 
shows that with a given incident angle and bare dielec- 



tric constants, we can deduce the information that Qf, is 



Fig. [4] (a) shows the ellipsometry angles comparing our 
approximated expression (Eq. ( |50[)) and full simulation 
results based on Eqs. ((33) and (1391). We can see that 



the linear approximation is valid even for high occupa- 
tion of xy and yz bands and also find that ignoring the 
contribution from ((AeTj, 1 )) gives results closed to the full 
simulation values. In Fig. [4] (b) the ellipsometry angles 
near w pe is investigated. As expected the reduction of the 
dip value is observed and additionally for high scattering 
rate (1/r ~ uj pe ) we find that a small peak emerges in 
the lower energy side of oj pe . 

We also investigate the peak visible in the higher fre- 
quency range on Fig. (a). In this case we can apply 

TT^, \AAV oi^^Q K^fV, n~Vr,„ 



Eq. (44), since both r^Tand r® give dominant contribu- 
tion to the reflectivity for the frequency away from ui pe . 
In this region, the bare reflectivity is mostly real and both 
it and the in-plane dielectric tensor (as seen in Fig.[3^b)) 
have smooth behavior. Thus we can expect that the peak 
in the ellipsometry angles are from the out-of-plane di- 
electric tensor ((Aej/)). To see this, we rewrite Eq. (44) 
up to the linear order of electronic dielectric tensor as 



24 r° 



(1 



(1 



„0\2 



M«Ae,.))] 



e° sin 2 6 



Im[(e&) a «Ae- 1 ))] 



(52) 



proportional to n 



• yz 



and the scattering rate 1/t. 



where we ignore the imaginary part of e^. Given that 
reflectivity are not sensitive to the frequency and both 
the sign of r° and r p are negative, we can see that the 
effect of in-plane dielectric tensor is to give smooth back- 
ground oc — f2 2 /(u; 2 + 1/r 2 ) and the peak is mostly due 
to the imaginary part of the out-of-the-plane dielectric 
tensor. The relation can also seen in Fig. [3] (b) where the 
peak energy of — ka.[qd{{Ae~g))] coincides with the peak 
in the right inset in Fig. [3] (a) for no = 0.3. 

In Sec. [n (b), we express the out-of-plane dielectric 
tensor with the density correlation function of xz/yz elec- 
trons (Eq. ( |23[ )). Since the imaginary part of the density 
correlation function represents the plasmon excitations, 
we fu rthe r investigate the density correlation function II 
(Eq. (22)). We can write the inverse of the II as 



„ 1 , , ui(oj + i/T)5 mm i — b^ nmo 8 mra i — 4n mo A. mmg v mm r A r , 

n mm'( w J - 7 / A fAT^ 

Emi U m mi + i l T ) - 9mA U mi m> 



(53) 
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where U mm > is a frequency independent unitary 
matrix satisfying U^ mi g mi U mim > = /\ 2 mmQ 5 mm , + 
4n mo y/ A mmo v mm / yj A m / mo using that the bare dielec- 
tric constant of STO does not change significantly at fre- 
quencies above the STO phonon frequency and is treated 
as a constant. Thus, we can get a simplified expression 
for ((Ae- 1 )) as 



( e ° s ) 2 «Ae z -V))) = ]T n } 1 ^ 

^— ' U}[U3 +1/T)- g r , 



(54) 



where f m = J2 m ' U mm >^2A. m > mo /t(m'\z\m Q ). 
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FIG. 5: (color online) (a) Diagonalized matrix elements (red 
dashed line) and density correlation function (blue line) for 
no = 0.3. (b) Comparison of ellipsometry angles with the 
different scattering rates for polar charge of 0.3 electrons per 
unit cell around plasmon peak, (c) Peak positions in ellip- 
sometry spectra compared with the peak in the out-of-plane 
dielectric tensor. 



With Eq. (54) we interpret the structure of the out- 
of-plane dielectric constant obtained by RPA. Without 
screening we have poles in II at energies corresponds to 
interband transition between the singly occupied state of 
xz (yz) band and the unoccupied states. The matrix el- 
ements are monotonically decreasing as increasing inter- 
band energy. However including the density response 



(screening) of the electron, expressed as Coulomb ma- 
trix Vram' , shifts the poles and mixes the matrix elements 
where shifted pole frequencies represent plasmon excita- 
tions related to the longitudinal oscillation of the xz / yz 
electrons and their oscillator strength arc determined by 
the matrix elements f m which are linear combinations 
of dipole matrix elements. Thus, we can interpret the 
peak of the ellipsometry angle as arising from plasmon 
excitations that are determined by two factors: how the 
eigenvalues g m are close to each other and the oscillator 
strength f m which are determined by band structure of 
the yz/xz electrons. In Fig. [6] (a) shows the diagonalized 
density correlation functions and its oscillator strength. 
We can see that the dominant contribution is from the os- 
cillator strength while the eigenvalue effects are relatively 
small. The rapid decrease of \f m \ 2 for y/g~^ > 0.25ey 
gives the asymmetric peak shape. The relation between 
plasma poles and the scattering rate is also investigated 
in Fig. [6] (b). For 1/t -C E Peak defined as the energy of 
plasma peak (Fig.[6](c)) the plasma poles are well defined 
and the scattering rate can be seen as a broadening in 
the plasma poles and it is expected that the peak height 
scales as f^r/u;. For 1/t 3> E Peak , the plasma poles no 
longer exist and the plasma peak vanishes. The plasma 
peak estimated by Eq. ( 54 ) is compared with ellipsom- 
etry data in Fig. [6] (c) where for high occupation, they 
agrees with each other and the difference is increases for 
smaller occupation since in Eq. ( 54 ) , we use dielectric 
constant of STO in the limit of u> — > oo. As the plas- 
mon peak energy decreases, the bare dielectric constant 
of STO deviates significantly from its limiting value and 
near oj pe the approximation used in Eq. ( 44 ) is not valid 
any more. 

We compare the position of the plasmon peaks with a 
simple estimation of the free electron plasma frequency 
fi£ = / 87re t(nxz+n V z) assoc i a t e d with an electron gas 

of volume density corresponding to our total areal den- 
sity of interface xz/yz electrons spread over a length L* 
which we take to be twice the peak position of the xz/yz 
density profile. The estimated plasma frequency for the 
polar charge density of {0.1,0.3,0.5} are {79,237,512} 
meV which are fairly close to the calculated plasmon 
peak position. The reasonable agreement between cal- 
culation and free electron model is expected since the 
binding energy of the xz/yz electrons is about two order 
of magnitude smaller than the plasma frequency. 

In Sec. [Tl] (c) presents detailed formulas for the re- 
flectivity by taking the long-wavelength approximation 
to the non-local dielectric tensor. The final results can 
be understood as implying that the interface system is 
described by an anisotropic dielectric tensor with com- 
ponents 



2ikd{{Ae {l (u))) 
^4) 2 <<Ae-V)>) 



(55) 



The complete expression of the effective dielectric tensor 
and its characteristic features are presented in Appendix 
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B. We note that the dielectric tensor depends on the in- 
cident angle through kd. The effective dielectric tensor 
allows us to have an alternative description of the inter- 
face electrons system. 

In this section, we investigated the features in the el- 
lipsometry measurements. By taking simplifying limits 
we are able to find the physical origin of the features and 
relate them with the physical quantities of interface elec- 
trons. Moreover we find the effective expression in terms 
of local dielectric tensor that may be helpful to provide 
another way to understand the behavior of interface elec- 
trons. Although the mapping between the plasmon peak 
and occupation of xz/yz elections and its scattering rates 
depends on the detail of band structure, combined with 
the information of the dip feature our results may allow 
us to reasonably estimate the physical properties of the 
interface electrons. 



V. SUMMARY 

We have presented an analysis of the subband struc- 
ture and optical response of an electron gas confined at a 
polar interface. We assumed that host material is similar 
to SrTiOa in having a conduction band derived from t2 g 
symmetry d-orbitals and possessing a very anisotropic 
dispersion and that the interface could be described as a 
charged sheet with areal charge density n which is a pa- 
rameter of the model. We used a self-consistent Hartree 
approximation to treat the Coulomb confinement and 
computed the dielectric tensor within the RPA approxi- 
mation. Our model involves several approximations. The 
most crucial of these is the neglect of mixing among the 
three t2g-derived bands. Actual band-mixing effects arise 
from spin orbit coupling and further-neighbor hopping 
and are small in the situation we consider 44 . In partic- 
ular, as discussed at length above, the response to in- 
plane electric fields is determined mainly by the number 
of carriers free to move parallel to the plane of the inter- 
face while the plasmon frequencies have to do with the 
total density of carriers with significant out-of-plane dis- 
persion; neither of these is affected much by the small 
interband couplings we have neglected. The neglect of 
band-mixing crucially simplifies the calculations, while 
the flexibility of our model-system formalism means that 
a wide range of situations can be studied and the sen- 
sitivity of the results to effects such as the total polar 
charge (controlled in part by defects which may be hard 
to model ab-initio) leads to increased insight. 

One crucial finding (implicit in previous work but ap- 
parently not remarked upon explicitly in the literature) 
is that the majority of the induced charge resides in xy- 
derived bands confined very closely (within 2 unit cells) 
of the interface. The charge density at farther distances 
has two components, with somewhat different distance 
dependences. One component arises from the lowest 
xz/yz subband, which is in most circumstances the only 
xz/yz subband to be occupied. The spatial decay of the 



charge density in this subband is determined by the sub- 
band wave function. The second component arises from 
a "Thomas Fermi tail" of electrons in xy-derived bands, 
with occupancy controlled by the self-consistent screen- 
ing of the interface potential. The charge density in this 
Thomas-Fermi tail is very small, but because of the very 
large dielectric constant of SrTi03 extends a very long 
distance from the interface. 

By use of our calculated subband structure, the RPA 
approximation and a Greens function formalism for the 
light propagation in an inhomogeneous situation we cal- 
culated the optical response of the interface system. Fol- 
lowing the experimental literature^ we presented the re- 
sults in term of changes in ellipsometry angles between 
the system with and without the interface electron gas. 
We find that the frequency dependent change in the ellip- 
sometry angle involves two important structures: a dip, 
at a frequency corresponding to the matching dielectric 
constants of SrTiOa and LaA103 and a peak at a higher 
frequency. The frequency at which the dip occurs is con- 
trolled by details of the bulk STO and LAO optical re- 
sponses and is not interesting. The magnitude of the dip 
was shown to be controlled by the in-plane conductivity. 
The peak in higher energy is in essence a c-axis plas- 
mon excitation of the interface electrons. We presented 
simple estimates allowing extraction of these parameters 
from the data without performing an elaborate fit of our 
calculations. 

Our results show how ellipsometry experiments can be 
a generally useful tool in the analysis of the electronic 
properties of oxide interfaces. Given the angle of the 
incident light and the bare dielectric constants of STO 
and LAO, our results shows that the location of the plas- 
mon peak and its height are easily mapped to the den- 
sity n yz and the relaxation time r through Eqs. (40) and 



( 52 ) . Moreover with the magnitude of dip at ui pe , and 



obtained r we can also obtain n xy + n yz (Eq. ( 50)). 
Thus the information about the orbital occupancy and 
disorders can be inferred. It has been reported the dis- 
crepancy between the number of carriers inferred from 
transport measurments^ and the nu mber expected from 
the polar catastrophe scenario ^ 26 * 38 ! or measured in high 
energy photoemission experiments^ and the mechanism 
that cause the discrepancy is not clear. Our results may 
provide insight to this by the comparison with the recent 
ellipsometry measurements^ that shows consistent dip 
and peak features. 
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Appendix A: Electromagnetic Green Function 

We introduce the electromagnetic Green function that 
is a solution of Eq. ( 27 1 following the derivation by Bagchi 
et a/.^. For s— and p— polarization, the explicit form of 
the Green function satisfying proper boundary conditions 
is presented. 

For a s-polarized light, the electric field is perpendic- 
ular to the plane of incident which is assumed to be in 
the xz plane. Thus, the electric field is only in the y 
direction. Assuming E oc ye l( ^ x E y (z) we can write 



V x (V x E) 



= V 2 E-V(V-E) 

= y(-Q 2 + ^)E y ( z y^ .(Ai) 



Then Eq. ( 25 1 becomes 



dz 2 



dziAe(z, zi)Ey(zi) 



(A2) 



q = 



1 _ 1 _ q + k 

2ik(l + rf) ~ 2iq(l - rf) ~ Aiqk 
q — k 



k 



€, LO COS ( 



k = 



sin 9 



(A7) 
(A8) 
(A9) 



and show that it satisfies Eq. ( A3 ) and proper boundary 
conditions. 

We can see that each of U y and V y satisfies the ho- 



mogeneous wave equation (Eq. (|28|). Therefore from 
Eq. (26), we expect that in the limit of z — > — oo, the 
asymptotic electric field becomes E y (z) cx e~ ikz , which is 
consistent with the expression of the self-consistent equa- 
tion in Eq. (29). Moreover, the continuity in U y and V y 
across the z — ensures that the E y is continuous across 
the interface. The interpretation of the Green function 
follows. For example, in the case of z, z' > the Green 
function is 



G yy (-^? Z ) 



lik \ 



\z-z'\ 



r O e ik(z+z') 



(A10) 



and the equation for the Green function is 



dz 2 



Q 2 G yy (z, zi) + —e°(z)G yy (z, z±) = 5(z - Zi) . 
We now claim that the Green function is following form: 

Gyy(z, Z) 



where 



U V (Z) = 



p[e(z - z')U y {z)V y {z') 

+Q( Z ' - Z)Vy(z)U y (z')] 



e iqz + r O e -tqz z < Q 

(l + r° s y kz z>0 



(A4) 



(A5) 



Viz) -I 0--r°)e-^' z<0 

V V\ Z > — 1 e -ikz _ r O e ikz z > ' ^ > 



The first term corresponds to the propagation of the wave 
from z' to z and the second term is the wave reflected at 
the interface and propagated back to the position z. For 



< and z 1 > 0, the Green function is written as 



Gyy(z,z') = —e-^(l + r° s )e 



ikz' 



(All) 



In this case, the Green function the wave propagate from 
the interface which is proportional to the left-moving 
wave from z' and previously reached wave which now 
reflected back to z' . 

Similar to the argument for s-polarized light, we can 
write differential equations for the Green function in p- 
polarized light case where electric field is in xz plane. 
The differential equation for the Green function becomes 
a coupled equation written as 



dz 1 



^e°(z) 



G xx (z,z) G xz (z,z') 
Gzx{ z , z ') G zz (z 7 z') 



S(z-z') 
S(z- 



(A12) 



The solution of the coupled equation is given as 



where 



G l3 {z,z>) = 



aj[e{z - z')U i {z)V j (z') 
+G(z' - zWzpjiz')] 



+ 6n5, 



•r iz e o 



e°(z)w : 



a { x z} 



-5(z - z') , (A13) 



Ai lj 2 / c 2 
1 q/* L +k/e° s 
Ai uj 2 / c 2 



c 2 q/e° L 
2iu) 2 l+r° 

c 2 g/ej 
2iu 2 l+r" 



fc/4 



k/e% 



(A14) 



16 



U x (z) 
V x (z) ■■ 
U z (z) 



e iqz _ r O £ -iqz z < q 

(1 - r°)e lkz z > 
(l + r°)e~ iqz z<0 

e -ihz + r O e ikz Z> Q 

iQ dU x (z) 
e°(z)ui 2 /c 2 - Q 2 dz 
Q [ q(e lqz + r°e-' l i z ) z < 

z > 



(A15) 
(A16) 



c°{z)uj 2 /c 2 -Q 2 \ fc(l-r£)e 



(A17) 



V z (z) 



iQ 



dV x {z) 



e°(z)cu 2 /c 2 -Q 2 dz 

q[l + r° p )e- iqz z<0 



(z)uj 2 /c 2 ~Q 2 \ k(e- tkz - ry kz ) z>0 



e° S q-e L k 
e s q + e° L k 



(A18) 
(A19) 



As in the case of s-polarized light, each set of {U X ,U Z } 
and {V^jVz} is the solution of Eq. (28), the homoge- 
neous wave equation. Thus, the Green function also sat- 
isfies the asymptotic boundary condition of electric field 
(E x oc V x ~ e~ iqz , E z oc V z - (Q/q)e~ lqz ) in the limit of 
z — > — oo. In Eq. (A13), there is one addition term for 



G zz components, which is required to give zero for the 
equation 



92 

dz- 



L0 2 €°(Z) 



G xz (z,z')-iQ-^-G zz (z,z')=0 



(A20) 

and Similar to the G yy in s-polarized light, we can write 
the Green function G xx for p-polarization case for z, z' > 
as 
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FIG. 6: (color online) Effective anisotropic dielectric constant 
for no = 0.3. (a) The real part and (b) imaginary part of 
anisotropic dielectric tensor compared with bare STO dielec- 
tric constant. Inset: magnification of the low frequency re- 
gion. 



Appendix B: Effective anisotropic dielectric tensor 



G xx (-3. ~ ) 



c 2 k 
2iu 2 e% 



e ik\z-z'\ _|_ r e ik(z+z') 



and for z > 0, z < as 



2 
c z q 



0\„ikz' 



(A21) 



(A22) 



which can be interpreted as the same way as the s- 
polarization case. For G zz (z — z') for z, z' > 0, we have 



G zz (z,z') 



Q 2 c 2 
2ik w 2 e° s 



e ik\z-z'\ _ r e ik(z+z') 



€i;UJ 



(A23) 



The delta function emerging in G zz accounts for the dis- 
continuity in the electric field due to the current in the 
z direction. We can eliminate the discontinuity by intro- 
ducing the displacement D z that is continuous across the 
interface. 



We introduce an anisotropic dielectric tensor en and 
e zz as an alternative description of the optical responses. 
The formulae are obtained by comparing the obtained ex- 
pression of reflectivity (Eqs. (331 and (39l) to the general 
expression of reflectivity in anisotropic media. 

The reflectivity of s- and p- polarized light for 
anisotropic materials with cubic symmetry are given byS^l 



en — e°r sin 2 1 



'e" T cose — \ e 



e° L sin 2 1 



cos 0/yft[ - yJ(i zz -e° L sin 2 0)/(e zz e||) 



cos 6/y/?[ + J (e zz - e° L sin 2 6)/{e zz e\\) 



,(B1) 



where e||(w) and e zz (uj) are anisotropic dielectric tensors 
for fields parallel and perpendicular to the interface, re- 



spectively. By the comparison with Eqs. (33) and (39), 
we obtain 
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|(w) = e ^-^M(A e ||»-^((A £ ||)) 2 

ej| + ^^sin 2 0((Ae||)) 2 
= 1 - 2ifed(£°)-i { e o e -||((A^ 2 1 )) + «A C ||»> + ^sin 2 ^e° ie -||((A e ^))2 ' lB2) 



We note that the obtained dielectric constants depend 
on the incident angle 9. We can think of en and e zz as 
an effective dielectric tensor which generates the same 
reflectivity as is given by our full calculation. Up to the 
first order of ujd/c, Eqs. ( |B2[ ) take the simple form, 

ey ~ e° - 2ikd{{Ae {] )) 
e zz ~ e° s + 2 l kd(4f((Ae- z 1 )) . 

Fig.[6]shows the real and imaginary part of the effective 
dielectric tensor for no = 0.3 with incident angle 75° 



compared with bare dielectric constant of STO. We can 
see the significant change in both the real and imaginary 
part contributed from the out-of-plane dielectric tensor 
around plasmon peak (~ 0.25eV). The in-plane dielectric 
tensor mostly contributes to the imaginary part of the 
dielectric tensor. In addition to the overall increases in 
the imaginary part of in-plane dielectric constant there 
is a dip near the to pe that is due to the rapid change in 
kd around LO phonon frequency of STO. The obtained 
dielectric tensor provides an alternative way to interpret 
the reflectivity obtained in Sec. [n] (c) . 
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